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Mining the secretome of Delia antiqua ( Diptera: Anthomyiidae ) 


based on the transcriptomic and gene expression data and 


analysis of its potential roles in diapause development 
HAO You-Jin* , XU Yan-Ling, CHEN Bin 


( Chongqing Key Laboratory of Vector Insects, Institute of Entomology and Molecular Biology, 


Chongqing Normal University, Chongqing 401331, China) 


Abstract; [Aim] Insect diapause is a critical eco-physiological adaptation for seasonal survival, but little is 


known about the molecular mechanisms of summer and winter diapause of Delia antiqua. Here, our main goal is 


to identify candidate genes linked to pupal diapause of D. antiqua. [Methods] A straightforward in silico 


approach was performed to predict the secreted diapause-related proteins based on the RNA-seq and digital gene 


expression datasets. Amino acid sequence similarity and basic physiochemical characteristics analysis were 


conducted. Gene expression levels of 12 selected genes from RNA-seq data were further validated in non- 


diapause pupae, and different developmental stages of summer and winter diapause pupae by qRT-PCR. 


[ Results] A total of 38 putative secreted proteins were differentially regulated between the summer and winter 


diapause pupae of D. antiqua. Gene Ontology analysis revealed differentially expressed genes to be associated 


with odorant binding, chitin and lipid metabolic processes, innate immunity and development, etc. Among 


them, regulation of lipid mobilization and chitin-related metabolism was found to be important during diapause 


development. qRT-PCR analyses showed that the expression profiles of the 12 selected genes were in good 
agreement with results from RNA-seq ( R =0. 862, P =0.001). [Conclusion] This study establishes a robust 
pipeline for the discovery of diapause-related genes by secretome mining based on the high-throughput 


sequencing. 
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1 INTRODUCTION 


The onion fly, Delia antiqua, is an important 
pest and widely distributes in Europe, Asian, and 
North America (Hao et al., 2016) and can exploit 
summer or winter pupal diapause to survive the 
severe seasons. Both summer and winter diapause 
occurs when the head evagination is completed. It 
makes D. antiqua an ideal insect model to decipher 
and compare molecular mechanisms regulating 
summer diapause (SD) and winter diapause (WD). 
It has been known that diapause is characterized by a 
down-regulation in the expression of many genes as 
well as the up-regulation of others ( Denlinger, 
2002). Recent molecular studies have uncovered a 
variety of diapause-specific genes from various insect 
species (Qi et al., 2015, 2016; Hao et al., 2016; 
Meyers et al., 2016; Zhao et al., 2017). Yet, the 
molecular and physiological 


precise processes 
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involved in the induction, maintenance and 
termination of diapause remain largely unknown. 
Therefore, it will be interesting to see whether D. 
antiqua exhibits a similar array of genes during its 
diapause. 

During insect development, some secreted 
proteins are essential for cell communications within 
the organism ( Horabin et al., 2003) or modulating 
responses to environmental stresses ( Durand et al., 
2010). Others have enzymatic functions involved in 
food degradation ( Robich and Denlinger, 2005; 
Macrae, 2010), and body remodeling ( Jia et al., 
2014) during the life stage transition. In this study, 
putative secreted proteins were predicted based on 
the first transcriptome of D. antiqua and gene 
expression was analyzed through seven digital gene 
expression ( DGE ) datasets. Interestingly, some 
gene homologues potentially associated with odorant, 
lipid metabolic innate 


chitin and processes, 
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immunity, etc., were firstly identified and in the laboratory for successive generations and 
differentially expressed at least one diapause diapause pupae were generated as previously 


developmental phase. These findings serve as an 
important basis for understanding the molecular 
mechanisms which underlies the pupal diapause 


development of D. antiqua. 


2 MATERIALS AND METHODS 


2.1 Insect rearing 

Delia antiqua were reared on artificial diets at 
50% -70% relative humidity under a photoperiod of 
16L: 8D at 20 +0.5°C, and these insects were kept 
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“Step I: Transcriptome and digital gene 


expression (DGE) data process 


Total RNA from ND, SD and WD 






RNA-S 
Short clean reads SS 


(51497228) 


Assembled ESTs (de novo SOAP) Short reads 
Unigenes Contigs DGE analysis 
(29659) (69176) 


Peptide sequences (ESTScan) 
(25152) 


described (Fu et al., 2016). Pupae were sampled 
according to the date (in days) after the pupation 
(DO), washed with ddH,O, and then stored at 
- 80. 
2.2 In silicon secretome prediction 

Nucleotide sequences were retrieved from the 
transcriptome data of D. antiqua ( SRR916227 ) 
(Zhang et al., 2014). Approaches to predict and 
annotate excretory proteins were performed by the 
following three steps (Fig. 1). 
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Fig. 1 


Bioinformatic flowchart for the secretome prediction 
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Step HT. Diapause-related gene annotation 


5 
DTN 


of diapause-related genes in Delia antique based on transcriptome 


and digital gene expression datasets from non-diapause pupae (ND), summer diapause pupae (SD) 


and winter diapause pupae (WD) 


The in silico analysis was performed with high stringent selection criteria. “Total number of short reads sequenced by next generation sequencing technique 


to explore onion maggot diapause mechanism. PNN = neural network. ° HMM = hidden Markov models. “TMHMM = trans-membrane hidden Markov 


model. ° Identified proteins are those proteins with N-terminus Met in Blast analysis. ‘Pineers are those proteins without N-terminus Met. 


2.2.1 
(DGE ) analysis: Transcriptome data and digital 
gene expression libraries were generated for non- 
diapause (ND), initial phase of summer and winter 
diapause (SI and WI), maintenance phase of 
summer and winter diapause (SM and WM), 


Transcriptome and digital gene expression 


quiescence phase of summer and winter diapause 
(SQ and WQ). Sequencing and assembling, and 


annotation procedures were described in a previous 
study (Hao et al., 2016). To create a list of 
differentially expressed genes with high confidence 
for further analysis, a stringent false discovery rate 
(FDR) (FDRS0.001) cutoff and | logg "=! | = 
1 were used as the threshold. 


2. 2. 2 


amino 


Excretory proteins prediction: Firstly, 


acid sequences retrieved from the 


1170 BR Acta Entomologica Sinica 60 ZS 





transcriptomic data of D. antiqua were fed into 
Signal P (http; // www. cbs. dtu. dk/services/ 
SignalP/) to predict the presence and cleavage site 
of signal peptides. Sequences with N-terminal signal 
peptide were passed to Wolf PSORT ( https: Ai 
wolfpsort. hec. jp/) and TargetP (http; // www. cbs. 
dtu. dk/services/TargetP/) for subcellular location 
prediction. The obtained dataset was then analyzed 
by TMHMM (http: // www. cbs. dtu. dk/services/ 
TMHMM/). To define a positive hit, the following 
criteria were used: (a) signal peptide predicted by 
neural network (NN) with the scores mean S and D; 
(b) signal peptide predicted by Hidden Markov 
models (HMM) considering the value of probability ; 
(c) Wolf PSORT (score cut-off =30) and TargetP 
(score cut-off = 95% ); and (d) signal peptide 
cleavage site located 10 — 40 amino acids from the N- 
terminus. Furthermore, because TMHMM may not 
distinguish signal peptides from transmembrane 
deduced 


transmembrane domain with 70 amino acids of the N- 


domains, proteins with single 
terminus were also considered as potential secreted 
proteins. 

2.2.3 Diapause-related gene identification and 
annotation; A unigene was considered to be a 
diapause-related gene when it produces a P <0. 05 


using the Chi-square algorithm in DGE analysis. 


Amino acid sequence similarities with known 
secreted proteins were searched using BlastP. The 
candidate sequences without transmembrane 


segments were passed to further analysis. Protein 


domain and family classification were analyzed by 
SMART (http; // smart. embl-heidelberg. de/) and 
Pfam (http; //pfam. sanger. ac. uk/). 
2.3 Digital gene expression validation by qRT- 
PCR 

To validate gene expression patterns identified 
in RNA-seq analyses, qRT-PCR was performed for 
12 selected genes in different function categories. 
Total RNA was extracted from one pupa/each 
sampling point with Trizol reagent ( Invitrogen, 
USA) according to the manufacturer’ s instructions. 
One pg total RNA was reverse-transcribed to 
generate the first-strand cDNA using SuperScript II 
reverse transcriptase (Invitrogen, USA). Each PCR 
reaction (20 uL in total) including 2 uL cDNA, 10 
pL 2x SYBR® Premix Ex Taq™ [I (Invitrogen, 
USA), 0.8 uL each of forward and reverse primer 
(10 pmol/L) was performed with the following 
conditions: 94°C for 2 min, followed by 40 cycles of 
95°C for 15 s, 55 — 58°C (dependent on gene of 
interest) for 15 s and 72°C for 20 s. The specificity 
of the amplicons were verified by melting curve and 
the fold change was calculated using the 27^“ 
method ( Livak and Schmittgen, 2001). Three 
biological replicates and three technical replicates 
were performed for each sample. The house-keeping 
gene GADPH was used as the internal control. The 
primers were designed based on the assembled 
nucleotide sequences in transcriptome data using the 
online server (http; // www. genscript. com/) and 
listed in Table 1. 


Table 1 The primers used for qRT-PCR 


Genes Forward primer (5' -3') 


Reverse primer (5' —3’) Amplicon size (bp) 





GCCTTTGTGAGGCAATTAGTGG 
CGTTGGGCCCGAATTGTGAT 
AAGGGCATGATTTGTGGTTCG 
ATGAATACTTGCCCGGCTGT 
CTTCCGCGATGACAGCTCTT 
TGAAATTGGGAGCACCGGCT 
ATTCTGCTGGTGCTTCACAC 
GCATGGTCGCCATAAAAGAGC 
TAATACGGGCAAGCCGTTCC 
GCAAAGGAGTATGTCGTGCC 
TCCCATGGTAGGTGGCCAAT 
AGCATTGCAACTTATGGCGT 
CGCCTTCTAGATACCGCCTG 


Destabilase 
Trypsin inhibitor 
Kazal type inhibitor 
Multiple domain protein 
Major royal jelly protein 
Chitinase 
Hypothetical protein 
Insect defensin A/B 
Cathepsin L 
Kunitz type inhibitor 
Cathepsin D 
Glycosyl hydrolase 
TAG lipase 


A RESULTS 


3.1 Assembly and extraction of transcriptome 
data 


Initially , clean reads 


51 497 228 short 


ACCGCCATCAGACCAAAAGG 101 
ATCATCTCTGGCAAACCCGT 113 
TTAGCAGGCTCCCTCCTTCT 123 
GCAGCGTGGAAGTACAAACC 142 
GGATTCATGGGTCCTACGGC 138 
TGGACCTTGGAAGGCAGTACC 122 
TTGGCGTTATTGCCGTTTGG 109 
ATTTCGGCACACGCAAACAG 139 
AACAGGAACCGCACACAGAT 125 
ATTTTCGTTGCCATCACAGCC 101 
CCATTTGGGCGATACGCAAG 131 
TCATTCCAATGGCCAGGACC 130 
AAATGCCAGCCGGATGAGTT 145 


(4 634 750 520 bases) with GC content of 39.03% 
were generated by deep RNA-seq from seven samples 
(ND, SI, SM, SQ, WI, WM and WQ) of D. 
antiqua. After SOAP de novo assembly, a total of 
69 176 contigs and 29 659 unigenes were obtained, 
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conceptually translated into 25 152 ploy-peptide 
sequences by ESTScan and used for the subsequent 
analysis. Additionally, 7 digital gene expression 
datasets were produced and mapped on the reference 
transcriptome. 
3. 2 Secretome and diapause-related gene 
prediction 

Firstly, 25 152 ploy-peptide sequences were 
analyzed by Signal P and generated 555 potentially 
secreted proteins. Then, 230 and 163 proteins were 
predicted as secreted proteins by TargetP and Wolf 
PSORT, respectively. A shared 150 sequences were 
passed to TMHMM for the 
transmembrane proteins. Thirty four proteins with 


prediction of 


one or more transmembrane fragments were removed 
from the secretory protein dataset. Ninety nine out of 
116 proteins with N-terminal methionine were 
Combined with 
digital gene expression data, finally 38 genes were 
predicted (Table 2). 
3.3 Sequence characteristics 

Five binding proteins including odorant binding 


predicted as secreted proteins. 


proteins (OBPs) , chemosensory protein (CSP) and 
chitin-binding proteins (Table 2) were identified, 
and sequence comparison showed that Ungiene 757 
and Unigene 2347 had 73% and 65% identity to 
OBP of respectively. 
However, they shared very low identity between each 


Stomoxys calceurans, 


other, implying they do not have a high conservation 
Unigene 674 
contains a conserved insect pheromone-binding 
domain ( Pfam; PF03392) and has highly conserved 
identity above 65% to Calliphora stygia CSPs 
( AID61324). 


Two insect 


and may have different functions. 


defensin A/B 
( Unigene103 and Unigene2131) were predicted and 
have 70% and 83% identity to their homologues of 
Protophormia terraenovae at the amino acid level, 


homologues 


respectively. However, only 34. 6% identity was 
shared between each other, suggesting that they do 
not display a high conservation. Domain analysis 
revealed that both have a typical insect defensin 
domain " knottins". Unigene780 encoding a lectin 
was significantly up-regulated in SM vs SI. Sequence 
analysis revealed that its full length amino acid 
sequence matched to C-type lectin of Drosophila 
willistoni (90% identity). Domain analysis showed 
that it contains a conserved carbohydrate- recognition 
domain (InterPro; IPR001304). 

Alignment analysis revealed that Unigene201 
shared 60% identity to a Kazal-type protease 
inhibitor of Drosophila erecta (EDV58866). While, 
Unigene828 appeared to be with full opening reading 


frame and shared 44% identity to a Kazal-type 
protease inhibitor of Drosophila grimshawi 
( EDV99366 ). Pairwise alignment revealed that they 
shared 37% identity to each other. In addition, 
Unigenel037 was predicted to be a kunitz-type 
inhibitor and shared 61% identity to the best- 
( NP _ 608799 ) in D. 


melanogaster at the amino acid level. Unigene 832 


matched homologues 


encoding a protein containing a typical serpin domain 
(Pfam: Pf 00079) and the amino acid sequence had 
67% identity to serpins of D. melanogaster. 
Interestingly, a cysteine protease inhibitor 
Unigene597 was identified and showed 54% identity 
to sarcocystatin of Sarcophaga crassipalpis 
(P31727). A conserved cystatin domain QxVxG (x 
can be of several amino acids) and a potential target 
protease binding sit G were found in amino acid 
sequence. 

Cathepsin D had 89% identity to their best hits 
in Musca domestica ( ABM69086 ) , while cathepsin 
L shared 91% identity to its homologue in 
Sarcophaga peregrina ( BAA76272 ). Sequence 
analysis revealed that both cathepsin D and cathepsin 
L consist of a signal peptide, a propeptide inhibitor 
domain (129) and a mature protein sequence. 

Proteins with putative functions differing from 
the above categories were grouped as proteins with 
energy homeostasis and miscellaneous functions. 
Transcripts without similarity to known genes were 
grouped as hypothetical proteins. A transcript 
( Unigene 2646) encoding a major royal jelly protein 
was identified and it had 75% identity to a 
of D. ( EDW37044 ). 
Similarity searching revealed that the full-length of 


homologue persimilis 
Ungiene5734 encodes a triacylglycerol lipase (TAG 
lipase) with 64% identity to its homologue of D. 
willistoni ( EDW76570 ). Protein domain analysis 
found an oa/B-hydrolase lipase domain ( Pfam; 
PF04083) and an a/f-hydrolase domain ( Pfam; 
PF00561). Unigene1435 encodes a multiple-domain 
protein including peritrophin A domain, low density 
lipoprotein receptor domain ( LDLR ), and chitin 
deacetylase-like domain. Similarity searching 
revealed that it has 92% identity to its homologues in 
D. mojavensis (EDW18929 ). 
3.4 Gene expression pattern analysis 

To gain insights into the global transcriptional 
patterns during the pupal diapause process, DGEs 
were identified from the normalized data by pairwise 
comparison of initiation phase of summer diapause vs 
non-diapause (SI vs ND) or initiation phase of 
winter diapause vs non-diapause ( WI vs ND), 


maintenance phase vs initiation phase (SI vs SM, 
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WM vs WI), and quiescence phase vs maintenance expressed at least one diapause phase and showed a 
phase (SQ vs SM, WQ vs WM). The fold changes very similar expression trend as in RNA-seq analysis 
were shown in Table 2. qRT-PCR analysis revealed with a spearman correlation of R = 0. 862 (P = 
showed that the 12 selected genes were differentially 0.001) (Fig. 2). 
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Fig. 2 Validation of RNA-seq results by qRT-PCR 
Expression levels of 12 selected genes in seven samples (ND, SI, SM, SQ, WI, WM and WQ) were detected by qRT-PCR. ND; Non-diapause; SI; 
Initial phase of summer diapause; SM; Maintenance phase of summer diapause ; SQ; Quiescence phase of summer diapause; WI; Initial phase of winter 
diapause; WM; Maintenance phase of winter diapause; WQ: Quiescence phase of winter diapause. Relative expression levels were derived from three 
biological replicates in qRT-PCR using ND as the calibrator (2 ~44“ method) , and log values of normalized RPKM of RNA-seq results. 














103 HAO You-Jin et al.: Mining the secretome of Delia antiqua (Diptera; Anthomyiidae ) 1175 





The transcripts of 9 genes were significantly up- 
regulated in SI vs ND (Table 2). The highest 
changed expression level was 8. 6-fold for tweedle E, 
then 5. 7-fold for destabilase, triacylglycerol lipase 
and new hypothetical protein ( Unigene 844). The 
transcripts of aminopeptidase (6. 5-fold) and two 
hypothetical proteins (5. 3-fold for Unigene 3981, 
and 3-fold for Unigene 844) were significantly 
higher at SM than at SQ. Whereas 8 genes were up- 
regulated in WI vs ND (Table 2). The expression 
level of a hypothetical gene ( Ungiene 844) was 
increased 59. 2-fold in WI vs ND, followed by 
Unigene 3981 (19. 4-fold) , OBP ( Unigene 2347 ) 
(15. 2-fold) , cathepsin D and cathepsin L (7. 3- 
fold) , and major royal jelly protein (4. 4-fold). In 
up-regulated in 
maintenance phase compared with the quiescence 


contrast, only 8 genes were 


phase. The most abundant transcripts were Unigene 
844 (9. 2-fold) , defensin A/B (Unigene 2131, 8- 
fold) and tweedle E (4.6-fold). 


4 DISCUSSION 


Based on a straightforward in silico approach, 
we generated the first comprehensive overview of the 
secreted proteins from the diapause onion fly pupa. 
Thirty eight genes were identified and potential 
functions of key genes during the diapause 
development were discussed as below. 

Insect OBPs are small, water-soluble proteins 
expressed in both olfactory and gustatory sensilla, as 
well as in other specialized tissues, and are secreted 
into the extracellular fluid. Generally, they were 
thought to be soluble hydrophilic odorants and bound 
by the respective receptors ( Andronopoulou et al., 
2006, Mitaka et al., 2011). Besides contributing to 
the recognition of odorants in insect, they also 
function as carriers in other developmental and 
physiological processes ( Foret and Maleszka, 
2006). As to the specific function of CSP and OBPs 
in the onion fly diapause pupae, we can only 
speculate. Perhaps the most intriguing hypothesis is 
that they might somehow regulate the chemosensory 
response to pheromonal odors in a hormonally and 
perhaps circadian sensitive manner. These findings 
served an important basis for understanding the 
molecular mechanisms of pheromones during the 
pupal diapause. 

Previous study also showed that defensin was 
up-regulated during larval diapause in spruce 
budworm Choristoneura fumiferana ( Palli et al., 
2001), but not in diapausing fly D.  triauraria 
( Daibo et al., 2001). In this study, the 


transcriptional level of defensin A/B ( Unigene 2131) 


was 8-fold higher in WQ vs WM. In contrast, it was 
down-regulated 7-fold in WI vs ND, and 13-fold in 
WI vs WM. This distinct expression patterns suggest 
that it could help winter diapause pupa increase 
ability to withstand attaches by pathogens or paly a 
crosstalk function of diapause development. The 
Drosophila gene fw (C-lectin orthologue ) mutation 
could severe impair eye and bristle development 
( Leshko-Lindsay and Corces , 1997). In flesh fly S. 
peregrina, a C-lectin was found to play roles in 
imaginal disc differentiation ( Kawaguchi et al., 
1991 ). 


regenectin ) was 


Another insect C-lectin ( designated as 
expressed selectively during 
cockroach tissue regeneration and proposed to play a 
potential role in control of development ( Arai et al., 
1998 ). The down-regulated expression pattern 
indicated that D. antiqua C-lectin might function 
directly or indirectly in immune defense and play a 
role in the development during the diapause. 
Destabilase is a member of lysozyme superfamily and 
its activity leads to the dissolution of stabilized 
fibrin. 


destabilase was identified in diapause insect and 


For the first time, a transcript coding 
significantly expressed at WI (Fig. 2). Due to its 
uncharacterized biochemical activities, it is still 
unknown whether this homologue has the destabilase 
activity or not and what is the real function in the 
winter diapause initiation. 

Two cystatins named Scys-A and -B have been 
previously described in the flesh fly S. crassipalpis. 
pupal 
diapause (Goto et al., 2002). However, Scys-B was 


Scys-A was down-regulated throughout 
highly expressed in early diapause. In this context, 
authors hypothesized that Scys-B would play a role in 
halting development. Interestingly, a similar 
expression pattern of Scys-B was observed for 
Unigene597. Its 


increased at WI, suggesting that it would have 


transcript was significantly 
potential role in initiating the developmental arrest. 


Being no cystatins have been biochemically 
characterized in diapause insects, what possible role 
it may play in D. antiqua is unclear. Serine 
proteases are involved in a variety of physiological 
functions such as development, digestion and stress 
responses. Loss of these critical protease inhibitors 
results in death of the organism or death. The up- 
regulation of a serpin ( HiDS) of the horn fly, 
Haematobia irritans has been observed and was 
thought to play an important role in inhibiting the 
differentiation and destruction of imaginal tissues 
during diapause ( Kalischuk-Tymensen, 2001). On 
the contrary, a serpin was down-regulated at the 


maintenance phase of the Colorado potato beetle, 


1176 BR Acta Entomologica Sinica 60 4 





Leptinotarsa decemlineata (Yocum et al., 2009). 
Here, D. antiqua pupae ceases of feeding, the 
diverse expression patterns of serpins between 
diapause and non-diapause suggests that they play 
some roles except the inhibiting the activity of 
digestive proteases. To our knowledge, serpins in 
diapause insects have not been functionally studied. 
What possible role serpins may play is unclear. 
Therefore , further experiments should be designed to 
investigate its real function in the pupal diapause of 
D. antiqua. 

Chitinases are members of the O-glycoside 
hydrolase superfamily and have been detected in 
molting fluid and gut tissues. In this study, Unigene 
3614 showed a U-shaped expression pattern in 
summer diapause pupae. However, its expression 
kept approximately constant across the whole winter 
diapause period (Fig. 2). As a state of 
developmental arrest, the 


lower expression of 


chitinase in diapause pupae of D. antiqua is 
understandable; however, its higher expression level 
at the quiescence phase suggests that it is associated 
with the termination of summer diapause and the 
preparation for adult eclosion. 

Imaginal disc growth factors (IDGFs) are a 
group of V chitinase-like proteins and have been 
shown to regulate cell proliferation and remodeling in 
insect. In D. melanogaster, four insulin-dependent 
IDGFs were found to promote cell proliferations 
(Kawamura et al., 1999). It was also identified in 
H. armigera and found to be associated with the 
pupal diapause (Chen et al., 2010). The highest 
expression of D. antiqua IDGF in ND pupae suggests 
that it plays a similar role in cell proliferations of 
imaginal disc. The decreased expression in the whole 
summer diapause period, WI and WM, and the up- 
regulated expression at WQ proved the arrested 
development in diapause process and development 
recovery in the quiescence phase. 

In this study, excreted proteins were predicted 
and their gene expression patterns were compared 
between different developmental stages to know about 
the putative transcripts involved in summer diapause 
and winter diapause. Many putative excreted protein 
genes showed no obvious similarity or convergence in 
two diapause types. More differentially expressed 
genes were found in maintenance phase of winter 
diapause than in summer diapause. Regulation of 
lipid mobilization and chitin-related metabolism at 
the transcriptional level was found to be important 
during diapause development and seasonal variation 
(summer or winter). The differences in the types 
and quantities of putative genes suggest that the 


evolution of two diapause types has led to the 
diversity of gene function adapting to specific 
seasonal environments. Taken together, our results 
provide a start point for the understanding of how 


excreted proteins function in the  diapause 


development of onion fly. 
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